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ABSTRACT 

We perform a timing analysis on RXTE data of the accreting millisecond pulsar 
XTE J 175 1-305 observed during the April 2002 outburst. After having corrected 
for Doppler effects on the pulse phases due to the orbital motion of the source, we 
performed a timing analysis on the phase delays, which gives, for the first time for this 
source, an estimate of the average spin frequency derivative < v >= (3.7±1.0) x 10~ 13 
Hz/s. We discuss the torque resulting from the spin- up of the neutron star deriving a 
dynamical estimate of the mass accretion rate and comparing it with the one obtained 
from X-ray flux. Constraints on the distance to the source are discussed, leading to a 
lower limit of ~ 6.7 kpc. 

Key words: stars: neutron - stars: magnetic fields - pulsars: general - pulsars: 
individual: XTE J 175 1-305- X-ray: binaries 



1 INTRODUCTION 

Accreting millisecond pulsars (AMSP in the following) are 
the long sought connection between low mass X-ray binaries 
(LMXBs) and millisecond radio pulsars. In fact, although it 
was hypothesised soon after their discovery that fast spin- 
ning radio pulsars were "recycled" by an accretion phase 
■ in a LMXB system, during which the neutron star (NS) is 
spun-up (see for a review Bhattacharya & van den Heuvel 
1991), evidence has been elusive since SAX J1808. 4-3658, 
the first accretion-driven millisecond X-ray pulsar, was dis- 
covered (Wijnands & van der Klis 1998). SAX J1808.4-3658, 
with a spin period of 2.5 ms, exhibiting both X-ray bursts 
and coherent pulsations, proved to be the missing link be- 
tween the two classes of sources. Since then, six more mil- 
lisecond X-ray pulsars were discovered (see Wijnands 2006 
for an observational review). 

All of these sources are transients with usually low duty 
cycles. Except for the case of HETE J1900. 1-2455 which re- 
mained active for more than a year after its discovery in 
June 2005 (Galloway et al. 2007), the outbursts of AMSP 
last for no more than a couple of months, with recurrence 
times usually larger than 2 yr (Galloway 2006). Although 
the sample is still small, monitoring of future outbursts ex- 
hibited by the known sources is extremely important for our 
understanding of LMXBs and their evolution. 
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The study of the rotational behaviour of these sources 
during outbursts is on the other hand obviously fundamen- 
tal as a test of theories of accretion physics. As a matter of 
fact, the measure of the variations of the spin frequency gives 
immediate understanding of the torques experienced by the 
compact object because of the accretion of matter, and fur- 
ther allows a model dependent dynamical estimate of the 
instantaneous mass accretion rate. Timing techniques per- 
formed on the coherently pulsed emission (see e.g. Blandford 
& Teukolsky 1976) represent the key tool in order to directly 
measure the variations of the spin rate of this kind of accre- 
tors. The application of this kind of analyses on the X-ray 
emission of these sources, as observed by the high temporal 
resolution satellite Rossi X-Ray Timing Explorer (RXTE) 
(Bradt et al. 1983), allowed the measurement of the spin fre- 
quency derivative in the case of IGR J00291+5934 (Falanga 
et al. 2005; Burderi et al. 2007), XTE J0929-314 (Galloway 
et al. 2002), XTE J1814-338 (Papitto et al. 2007), XTE 
J1807-294 (Riggio et al. 2007) and of one of the oubursts 
shown by SAX J1808.4-3658 (Burderi et al. 2006). See Di 
Salvo et al. (2007) and references therein for a review. 

Even though the shortness of the outbursts generally ex- 
hibited by AMSP strongly limits the capability of the timing 
analysis in discriminating between various accretion models, 
observations have already shown how the behaviour of these 
sources can be variable, in fact, as a result of accretion, some 
of them are observed to spin up, while others decelerate. In 
this work we present a timing analysis performed on the only 
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outburst of XTE J 175 1-305 observed so far by RXTE with 
its highest temporal resolution mode, and discuss the ob- 
served spin frequency evolution as a result of the accretion 
torques acting on the NS. Moreover the measure of the spin 
frequency derivative is used to infer a dynamical estimate 
of the peak mass accretion rate during the outburst, which 
can be compared with the estimate deduced from spectral 
modelling of the X-ray emission, in order to give an estimate 
of the distance to the source. 



2 OBSERVATIONS AND DATA ANALYSIS 

XTE J1751-305 was first detected monitoring the Galactic 
bulge region with the Proportional Counter Array (PCA) 
on board the RXTE (Markwardt & Swank 2002a). Subse- 
quently, pointed observations in April 2002, allowed the de- 
tection of an X-ray periodic modulation at a frequency of 
about 435 Hz, establishing that this source belongs to the 
class of accreting millisecond X-ray pulsars (Markwardt et 
al. 2002b, M02 hereafter). Throughout this paper we used a 
sample of RXTE public domain data taken between 4 and 
30 April 2002. The data we use to perform the timing anal- 
ysis are those collected by the Proportional Counter Ar- 
ray (PCA, Jahoda et al. 1996), which is made by five units 
(PCUs) sensitive in the 2 — 60 keV band, for an overall col- 
lecting area of ~ 6250 cm 2 . Except the observations of April 
4, which were taken in Good Xenon configuration with a 
time resolution (2 s read out time and 256 energy channels), 
all the other PCA data we considered are Generic Event 
(E_125us_64M_0_ls) with a temporal resolution of 122/is and 
64 energy channels. All the data were processed and anal- 
ysed using the HEASARC FTOOLS v.6.0. 

The source was first spotted by the PCA in the desired 
collecting mode at To = 52368.653 MJD, which represents 
the start time of observations throughout this paper. There- 
after the X-ray light curve of XTE J1751-305 is made of two 
exponential decays with different e-folding factors. The first 
one describes the first 8.5 d of observation and can be mod- 
elled with the function L x (t) = L x (T)exp{-(t - T)/t£ ] ], 
where L X (T) is the luminosity at f = 52369.644 MJD and 
t£ ] = 7.2 d (see Gierlinski & Poutanen (2005), GP here- 
inafter, for the 2 — 20 keV light curve; see also M02). Sub- 
sequently the light curve experiences a sharp break and can 
be described with a similar decay function this time with 
T< d^ ~ 0-^3 d. The source then switches back to the quies- 
cence emission levels, ~ 10 d after the first available observa- 
tion. Considering a physically motivated spectral model (see 
the discussion for further details) GP estimated for the bolo- 
metric X-ray/7-ray luminosity attained by XTE J1751-305 
during the outburst L x (f) = 2.7 x 10 37 d|. 6 erg s" 1 , where 
ds.5 is the distance to the source in units of 8.5 kpc. 

The technique we used in order to perform the timing 
analysis on the pulsed emission is extensively described in 
Burderi et al. (2007, see also Papitto et al. 2007). We first 
corrected X-ray photons arrival times to the Solar system 
barycentre considering the best position available for this 
source from Chandra (M02). We then focused on the orbital 
modulation of the phases to derive an accurate orbital solu- 
tion. The evolution of the phases, measured by folding 90 s 
long intervals around the M02 estimate of the spin period, 
was modelled with Eq.l of Papitto et al. (2007), without 



Table 1. Orbital and timing parameters of XTE J1751-305 



M02 



This work a 



Orbital solution 

a sin i/c (lt-ms) 

Porb (s) 

T* (MJD) 
Eccentricity e 



10.1134(83) 
2545.3414(38) 
52368.0128983(87) , 
< 1.7 x 10~ 3 



10.125(5) 
2545.342(2) 
52368.0129023(4) 
< 1.3 x 10" 3 



Timing Solution 

vo (Hz) 

< v > (Hz/s) 

!>(f ) (Hz/s) a = 2/7 



435.317993681(12) 
< 3 x 10~ 13 c 



435.31799357(4) 
(3.7 ± 1.0) x 10- 
(5.6 ± 1.2) x 10- 



a Numbers in parentheses, referring to our values, are the 90% 
confidence level uncertainties in the last significant figure, while 
the ones referring to M02 are given at 3<r confidence level on the 
last significant digits. The same confidence levels are considered 
in giving upper limits on the eccentricity, e. 

^ The value reported here is revised with respect to the one 
originally quoted in M02 (Markwardt et al. 2007). 
c This upper limit is to be considered on the absolute value of 
< i> >. 



considering the term owing to the position uncertainties (see 
below). The orbital parameters we obtain, namely the pro- 
jected semi-major axis measured in light ms, a sin i/c, the 
orbital period, P or b, the time of passage of the NS at the 
ascending node of the orbit , T* , and the eccentricity, e, are 
listed in TabfTJ which contains for comparison purposes also 
the ones reported by M02. The two orbital solutions are in 
good agreement whitin the quoted uncertainties. 

The times of arrival of X-ray photons were then re- 
ported to the line of nodes of the binary system orbit con- 
sidering our orbital solution. As the time over which any un- 
certainty on the orbital parameters may possibly affect the 
phases of the X-ray pulsations (P or b) is much smaller than 
the time required for the spin frequency derivative to pro- 
duce a significant effect, the two effects on residuals can be 
decoupled. The effect on pulse phases due to the remaining 
uncertainty on the orbital parameters, o> or (,, can therefore 
be treated as a normally distributed source of error (see Eq.3 
of Burderi et al. 2007 for an expression of <7^ or f>). The final 
uncertainty, on the phase residuals we consider in the 
following is then the squared sum between a$ or b and the 
statistical error arising from sinusoidal fitting of the pulse 
profiles. 

By folding each light curve corrected for the orbital mo- 
tion of the system around our guess for the spin frequency 
vf (which initially was the M02 estimate), we could detect 
coherent pulsations until MJD 52377.6, 9 d after the first 
observation available, when the X-ray flux had become ap- 
proximately one tenth of the peak flux. 

The criterion we considered in order to assess the pres- 
ence of a periodic signal is simply based on the expected 
statistic distribution followed by a folded light curve in 
which no signal is present (Leahy et al. 1983). The statistics 
5* used to check the presence of a periodic signal is defined 
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Figure 1. Plot of the observed statistics S b s of every consid- 
ered folded light curve. The horizontal line represents the 99% 
confidence limit for the detection of pulsations. 



< 

o 

Q_ 



430 



420 



410 - 




0.5 



1 

Pulse phase 



1.5 



Figure 2. Pulse profile produced folding in 20 phase bins around 
uq = 435.31799357 Hz, all the available observations spanning the 
time interval from MJD 52370.537 to MJD 52370.776 (OBSID 
70131-01-02-00). This time interval does not represents a contin- 
uous time series and the effective total integration time is ~ 14 ks. 
The solid line is the best fitting sinusoid. The reduced chi squared 
of this model is \ 2 = 0.78 and is not significantly improved by 
adding an harmonic. Two cycles are plotted for clarity. 



as 

s = Y. {Rj ~ 2 R) \ a) 

3=1 

where n is the number of phase bin by which the light curve 
is divided, j — 1, n, Rj is the count rate in the j-th bin, 
R is the average count rate in the considered observation, 
and a 2 = R x n/T, with T total time of integration. It 
is assumed that in the absence of periodicity the photon 
counting in each bin of a folded light curve, Rj x T/n, is pure 
Poissonian counting noise, so that its mean and variance can 
be estimated from the average number of counts in each bin, 
R x T/n. In the limit of a large number of counts, as in the 
case of the 2002 outburst of XTE J1751-305, when R x 
T/n Si 2000, we can assume that this distribution behaves 
as a Gaussian with the mean equal to the variance. If these 
conditions are met the statistic S of a folded light curve with 
no signal is expected to follow a chi-squared distribution 
with n — 1 degrees of freedom. It is then immediate to define 
a threshold for the detection of the signal So from the desired 
confidence level P, (1 — P/100) = Qn-i(Xo = "So), where the 
term on the right hand side is the integrated probability from 
Xo to co of a chi-squared distribution with n — 1 degrees of 
freedom. In this case we chose P = 99 and n = 19, so 
that every observation with Sobs < 36.19 was considered as 
containing no signal at the considered level of significance, 
and therefore withdrawn from the sample used in timing 
analysis (see Fig[TJ. 

Folded light curves that met the detection criterion were 
then modelled with a sinusoid with the period fixed to the 
folding one. We also tried to add an harmonic to the fitting 
function, but such a component was significantly detected 
and led to a slightly better quality fit only in a small fraction 
of the observations (^ 10%), so that we consider for the 
purposes of the timing analysis only the phases obtained 
with a single harmonic sinusoidal function. An example of 
the pulse profile is shown in Fig [2] 



Therefore we modelled phase residuals with a parabolic 
function 

A$(t) = 1/pf x{S + /3t + -/t 2 ), (2) 

where vp is the folding frequency, f3 — — Avo is the correc- 
tion to vf to obtain an estimate of the spin frequency at 
t = To, 7 = — < v > / 2 is the term owing to a constant 
spin frequency derivative and 8 is a phase constant. This 
procedure was repeated several times, correcting each time 
the frequency over which we fold the light curves, until the 
linear term of Eq[2] was compatible with being null within 
the uncertainties. The timing parameters we finally obtain 
are v = 435.31799357(3) Hz and < v >= 3.7(8) x 10" 13 
Hz/s, where the numbers in parentheses are the 90% confi- 
dence level uncertainties on the last significant figure, as for 
all the uncertainties quoted in the rest of the paper. 

In Fig[3]we plot the phase delays A$(£), measured in 
microseconds, of light curves folded around Vo, versus the 
time elapsed since the first observation considered. The plot- 
ted error bars refer to the overall la uncertainties on phases, 
computed by considering both the statistical errors coming 
from the sinusoidal modelling and the errors induced by the 
uncertainties on the orbital parameters listed in the right 
column of TabfTJ as already stated. The reduced chi squared 
of the quadratic best fitting model is 81/56, which compared 
to the one obtained with a linear fit (136/57), gives an F test 
probability of 8 x 10~ 8 that the improvement in the variance 
is purely given by chance. 

As already shown for other sources of this class (Burderi 
et al. 2006, Papitto et al. 2007, Riggio et al. 2007) the pres- 
ence of timing noise, probably due to changes in the position 
of the emitting hot spot on the NS surface, suggests cau- 
tion in considering the reliability of the uncertainties on the 
timing parameters obtained by a simple least square fit to 
the phase residuals, taken with their uncertainties a$. Nev- 
ertheless the pulse phase evolution of XTE J1751-305 can 
be modelled satisfactorily by a general parabolic trend (see 
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Figure 3. Evolution of the pulse phase delays, measured in 
microseconds, folding every available observation around vq = 
435.31799357 Hz. The plotted error bars are the la uncertainties 
on phases, including both the statistical errors from the pulse 
profile sinusoidal fitting and the errors induced by the residual 
uncertainties on orbital parameters. We note that these error bars 
have not been rescaled by the factor 1.2 that makes X? = 1> which 
we used to get the errors in the parameters given in Tab[TJ (see 
text for details). The solid line is the best fit constant i> model 
(upper panel). Residuals with respect to the best-fitting model 
(lower panel). 



Discussion for a comparison with the behaviour displayed 
by other sources of this class) , even if the presence of small 
timing residuals leads to a reduced \ 2 slightly higher than 
one. Therefore, in order to get a conservative estimate of the 
uncertainties affecting the measured spin frequency and its 
derivative, we amplified all the errors on our phase measure- 
ments by a common factor (1.2) until we obtain \\ ~ 1 from 
a fit with Eq[5] We then consider as the most conservative 
the uncertainties a VQ and so obtained. In this way we get 
to our final estimates of the timing parameters, which are 
the ones listed in Tab[TJ For sake of clarity we note however 
that the error bars plotted in Fig(3] are not multiplied by 
any factor and they represent the genuine la uncertainties 
on the measured phases. 

A more physically motivated description of the evolu- 
tion of phase delays is easily obtained if the dependency 
of v(t) on the accretion rate is considered. According to the 
simplest model of angular momentum exchange between the 
accreting matter and the compact object, matter in the disc 
can be considered fully diamagnetic, so that the only rel- 
evant torque is the positive one coming from accretion of 
matter at the inner boundary of the disc, which is placed 
at a distance Ri n from the centre of the compact object. In 
an X-ray pulsar this boundary can be determined, at least 
by orders of magnitude, considering the balance between 
the ram pressure of the infalling matter and the magnetic 
pressure of the field that truncates the disc. When mass 



accretion rate declines, as in the fading parts of the out- 
burst light curves of X-ray transients like XTE J1751-305, 
the pressure of falling matter is expected to decrease letting 
the magnetosphere to expand. Assuming that the inner disc 
boundary scales as Ri n oc M~ a , with a ^ and constant at 
every experienced accretion rate, the integral of the spin up 
torque exerted by matter accreting at Ri n can be therefore 
carried on in terms of the instantaneous value of the mass 
accretion rate. The instantaneous spin frequency derivative 
can then be expressed as (see e.g. Burderi et al. 2007): 



-±-l M(t) (W 



l-a/2 



(3) 



where lo = (GMRin(i)) 1 ^ 2 is the specific angular momen- 
tum of accreting matter evaluated at t = i . Once it is as- 
sumed that M(t) oc Lx(t), the temporal dependence of z>(t) 
can be expressed in terms of the particular shape of the out- 
burst light curve. As our aim is to compare our estimate of 
the accretion rate needed to explain the observed v assuming 
a certain torque model, and the accretion rate deduced from 
spectral fitting, we consider here and in the following i — T, 
an instant for which GP produce an accurate estimate of the 
0.7 - 200 keV flux. 

Considering an X-ray light curve shape composed only 
by a single exponential decay, the expected behaviour of the 
phase delays time evolution is then: 



A — [B + kr a C] (t - T ) - CtI exp 



t - T 



A$(t) = — 

(4) 

where r Q = r<j(l — a/2) , Td is the e- folding factor of the 
light curve exponential decay, k = exp [—(To — T) /T a ] , A is a 
phase constant, B = Avo is the correction to the folded spin 
frequency and C = u{T) = {2n I)- 1 [GMR ln (T)] 1/2 M(T) is 
the spin frequency derivative at t = T. 

We first considered eq.Q with Td = 7.2 d, and fitted the 
observed phase delays with various values of a, correspond- 
ing to different models of disc-magnetosphere interaction. In 
particular, the case a = corresponds to no dependence of 
Rin on the instantaneous accretion rate and in this case we 
consider Ri n = Rc where Rc is the corotation radius de- 
fined as the radius at which the Keplerian orbiting matter in 
the disc exactly corotates with the NS magnetosphere (see 
e.g. Rappaport et al. 2004). On the other hand, a — 2/7 is 
the Alfvenic value obtained in the approximation of spheri- 
cal accretion of matter. Unfortunately we succeded in having 
no significant improvement in the description of the phase 
residuals by implementing such torque models, as for exam- 
ple we obtain Xr = 79/56 when a = 2/7 is considered. 

We also tried other a plausible values, based on a more 
realistic treatment of the disc structure at the inner rim, 
which includes different possible regimes, as gas pressure 
dominated or radiation pressure dominated SS optically 
thick disc (respectively a = 0.25 and 0.15 for models 1G 
and 1R in Ghosh 1996 and Psaltis & Chakrabarty 1999). 
We did not consider any other model that would lead to 
significantly higher values of a (e.g. two temperature opti- 
cally thin gas pressure dominated discs), as the maximum 
value that this parameter can attain is the one that im- 
plies a full excursion of R in during the outburst, from the 
NS radius to the corotation radius, while the source keeps 
showing pulsations (see discussion). Since in the case of XTE 
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J1751-305 a coherent signal can be detected as long as the 
X-ray flux, which is related to the instantaneous mass accre- 
tion rate, has declined to one tenth of its maximum value, 
the constraint Ri n (tcut) / Rin(T ) = (F x (t cut )/F x (T ))- a < 
Rc/Rns, where t cut — To + 9 d is the time at which the pul- 
sations are cut off, leads to a < 0.4 for a 1.4 Mq NS with a 
Rns = 11.1 km radius. The available statistics nevertheless 
revealed itself to be too low to discriminate between these 
models. This is witnessed by the fact that the goodness of 
the fit and also the value of the spin frequency derivative 
remain almost the same when different plausible values of 
a are chosen. In order to give a reference value of u(T), the 
spin frequency derivative at t = T, we consider the simple 
Alfvenic case, a = 2/7, for which v(T) = (5.6 ± 1.0) x 10" 13 
Hz/s. In the same way we operated before, we amplified the 
actual errors on phases by a common factor until we have 
Xr — 1 for the best fitting a = 2/7 model, thus obtaining 
the 90% confidence level uncertainty on v(T) listed in Tab[T] 
We note that an attempt was also made to model the 
behaviour of the phase delays in terms of the broken decay 
effectively observed in the X-ray light curve, rather than 
considering a single decay taking place throughout the out- 
burst. We thus implemented in Eq(4] the change of the e- 
folding factor that takes place simultaneously to the break. 
This attempt anyway led neither to improvements in the 
quality of the fit nor to variations in the measured parame- 
ters for every value of a considered. This is probably due to 
the shortness of the time interval in which the light curve is 

(2) 

described by before the pulsations fade away. 

We conclude this section by noting that a systematic 
term due to the uncertainty in the source position has to be 
considered in order to have a reliable estimate of the uncer- 
tainty in the spin frequency as well as in the spin frequency 
derivative. The best position available has a 90% error circle 
of 0".6 (M02), and following the expression given by Bur- 
den et al. 2007, we estimate the upper limit on the effects of 
this uncertainty on the value of vo as cr sy st v () < 1 x 10 -7 Hz, 
while on i/ (T) as a syst c < 0.3 x 10" 13 Hz/s, which is one or- 
der of magnitude smaller than the spin frequency derivative 
we measure, and smaller than the error on < i> > quoted 
above. 



3 DISCUSSION 

In the previous section we described the application of dif- 
ferent accretion models, referring to different assumption on 
the disc-magnetosphere interaction and on the particular 
structure of the inner region of the accretion disc, to the 
observed evolution of the pulse phases of XTE J1751-305, 
under the assumption that the curvature observed in the 
phase delays versus time is indeed a measure of the spin 
frequency derivative. 

The stability of the pulse phase evolution of this source 
resembles the one shown by IGR J00291+5934 (Burderi et 
al. 2007) and XTE J0929-314 (Galloway et al. 2007) with 
a smooth variation along the course of the outburst, wit- 
nessed by the small post fit residuals obtained even when 
a simple constant v model is used. This behaviour is dif- 
ferent from the one analysed in an other subset of sources 
of this class, namely XTE J1814-338 (Papitto et al. 2007) 
and XTE J 1807-294 (Riggio et al. 2007), where the phases 



oscillate around the mean trend clearly anticorrelating with 
rapid (~ 1 d) X-ray flux variations, and SAX J1808.4-3658 
(Burderi et al. 2006), which shows an even more complex 
behaviour. It has to be noted that for all of the latter at 
least two harmonics are needed to model adequately the 
pulse profiles, while, among the sources of the first group, 
only XTE J0929-314 has a significant harmonic content. In 
particular the trend followed by the phases of the second har- 
monic seems more stable than the one of the fundamental 
(Burderi et al. 2006, Riggio et al. 2007), leading the authors 
to consider the second harmonic to establish the rotational 
behaviour of the considered source. A correlation between 
the phase stability and the two class of AMSPs depicted 
above seems to arise if it is considered that for the last three 
the X-ray light curve is somehow complex, showing a variety 
of features, like re-flares, oscillations around a mean trend 
and in general variations on time scales of days, and simul- 
taneously the pulse phases deviate from a continous trend. 
The light curves of the first three instead generally show 
a quite smooth exponential decay, during which the phases 
distribute themself normally around the mean trend. 

A discussion of the noise of AMSP pulse phases around 
their mean trend is beyond of the scope of this paper, never- 
theless these considerations make us even more confident on 
the stability of the evolution of the pulse phases in the case 
of XTE J1751-305, a source in which no secondary harmon- 
ics appear nor the light curve deviates from an exponential 
decay through the time interval considered here. We there- 
fore ascribe the parabolic trend followed by the phases to a 
spin frequency evolution due to the accretion torques acting 
on the NS when it is efficiently accreting mass. 

However the attempts made proved unsuccessful in dis- 
criminating a constant spin up model from a model in which 
the spin frequency derivative depends on the instantaneous 
value of the mass accretion rate, as traced by the observed 
X-ray flux. In the hypothesis that i>(t) oc M(t) tx Fx(t), 
which we stress is not favoured nor disfavoured by the data, 
the measurement of the spin frequency derivative of XTE 
J1751-305 at t = T, v(T), allows a dynamical estimate 
of the mass accretion rate. This can be further compared 
with the estimate of the flux observed jointly by XMM and 
RXTE, in order to constrain the distance to the source. The 
evaluation of Eq.(|3} at the time t — T, together with the 
estimate of is(T) given in the previous section, leads to an 
expression that relates mass accretion rate at that time to 
the lever arm of the spin up torque at the same time, i.e. 
the inner disc radius: 

M = (30 ± 6) x 10- 10 / 45 m- 2/3 <r 1/2 M Q /yr, (5) 

where J45 is the moment of inertia of the compact object 
in units of 10 45 g cm 2 , m is the mass of the NS in solar 
units, and £ = Rin/Rc is a parametrisation of the inner 
disc radius at T in terms of the corotation radius Rc- For 
XTE J1751-305 we have R c = (GM/Q%) 1/3 = 26.1m 1/3 
km, where M is the mass and ils is the angular rotational 
velocity of the neutron star. The use of such a parameter 
is particularly suited for an accreting pulsar, as it has to 
be Rns/Rc < < 1 in order for pulsations to be visi- 
ble at a certain time t. The lower constraint is an obvious 
consequence of the interpretation of coherent pulsations as 
due to funnelling to the magnetic poles of the accreted mat- 
ter, while the upper limit is due to the centrifugal inhibi- 
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tion of efficient accretion of matter onto the NS when the 
NS-magnetosphere system rotates faster than matter at the 
inner boundary of the disc. Expressing these boundaries in 
terms of m and Tie, the radius R of the NS in units of fO 6 
cm, one obtains for XTE J1751-305, 0.387? 6 m~ 1/3 < £ < 1, 
which allows the definition of a range of possible accretion 
rates at t = T, by their substitution in Eq. (|5jl : 

(30±6)/ 45 m- 2/3 < Mio(To) < (48±9)/ 45 7? 6 " 1/2 m- 1/2 (6) 

where M\q is the mass accretion rate in units of 10~ 10 
M /yr. 

This estimate of the mass accretion rate can be ex- 
pressed in terms of X-ray luminosity via the usual relation 
Lx = eGM M I Rns , where e ~ 1 is the rate of conversion of 
gravitational energy released in accretion to observable X- 
ray luminosity. This allows to make a comparison of the dy- 
namical estimate of the accretion luminosity we derived from 
timing analysis with the one obtained by spectral modelling 
of the observed X-ray flux, which is obviously dependent on 
the source distance. By considering the expression given by 
GP for the bolometric luminosity, Lx(T), one obtains for 
the distance d of XTE J1751-305: 

d = 8.2 I^m^R-^e^C'^kpc (7) 

Considering a moderately stiff EoS for an m = 1.4 NS, 
such as the FPS, one has for a compact object spinning 
at the rate measured for XTE J1751-305 7 45 = 1.24 and 
R& = 1.11 (see e.g. Cook, Shapiro & Teukolsky 1994). ft 
is then possible to find a range for the distance d of XTE 
J1751-305, 9.lkpc < de~ 1/2 < 11.6/cpc. 

XTE J1751-305 is located only 2° away from the Galac- 
tic Centre, so that a distance higher than 8.5 kpc would be 
highly improbable, if the source shows significant emission at 
energies of the order of ~ 1 keV (GP estimated Nh ~ 10 22 
cm -2 ). However, the discrepancy between our lower limit on 
the distance, 9.1 kpc, and the one implied by this constraint 
is not large, especially considering the relatively large uncer- 
tainties underlying the simple arguments which led to our 
distance determination, and the assumptions made on the 
torque model and on the NS structure. Moreover we recall 
that, because of the uncertainty in the source position, the 
estimate of u(T) we used to get a dynamical estimate of the 
mass accretion rate is affected by a systematic error, which 
can in principle decrease our estimate of the distance below 
8.5 kpc. M02 estimated a lower limit on the distance of 7 kpc 
based on indirect arguments, and our comparison between 
the observed v and the measured X-ray flux strongly sup- 
ports the hypothesis that the source is close to the Galactic 
Centre. 

Nevertheless, we can argue that the limits on the dis- 
tance we find above may be overestimated due to several 
reasons, among which, a non isotropic emission from the 
source or an occultation of part of the accretion luminos- 
ity. In particular, as already noticed by GP, the fact that 
the pulse shape shown by XTE J1751-305 is almost sinu- 
soidal implies that only one spot is visible from our line of 
sight. If this had not been the case, i.e. if also an antipo- 
dal spot had intercepted our line of sight during rotation, 
the folded light curves should have shown a secondary max- 
imum (or a plateau in a more extreme case), that would 
have required the addition of at least an harmonic to be 
efficiently modelled. This picture is consistent with the re- 



sults of detailed 3D magnetohydrodynamics simulations of 
disc accretion to this kind of rotators performed by Kulka- 
rni & Romanova (2005). We argue that mass is accreting on 
both the polar caps of the compact object, with a fraction of 
the emitted accretion luminosity being blocked by an opaque 
absorber (probably by the accretion disc, as the absence of a 
significant Compton reflection component strongly suggests 
i > 60°, see GP), and re-emitted outside the considered 
0.7 - 200 keV energy band (see also Burderi et al. 2007). It 
could be therefore the case that the value given by GP for 
the X-ray flux represents an underestimate of the real emis- 
sion owing to the accretion of matter. This in turn would 
decrease our distance estimates to values in a better agree- 
ment with the constraint of the source not being farther 
than the Galactic Centre. Defining rj as the ratio between 
the effective accretion luminosity and the one observed, we 
find lower distance estimates that would place the source not 
farther than the Galactic Centre for r\ — 1.2, and a lower 
limit of 6.7 kpc corresponding to r\ = 2, the maximum value 
r\ may reasonably attain. 

It has to be highlighted that all these estimates of the 
source distance are derived by assuming that no negative 
torque is acting on the NS, due e.g. to threading of the disc 
by the magnetic field lines in regions where matter in the 
disc rotates slower than the threading field lines, as it is as- 
sumed to explain the rotational behaviour of other source of 
this class (see e.g. Papitto et al. 2007) . If this would not have 
been the case, higher estimates of the mass accretion rate 
would be obtained, as the observed spin frequency deriva- 
tive would represent the balance between a negative torque 
and the positive one due to accretion at Ri n . This would 
straightforwardly lead to higher estimates of the distance. 



4 CONCLUSIONS 

We performed a detailed timing analysis upon the coherent 
pulsations shown by XTE J 175 1-305 during its 2002 out- 
burst as observed by RXTE. We could detect such a signal 
in the period 2002 April 4 — 13, i.e. in all observations during 
which the source emitted an X-ray flux above the quiescent 
level. 

After having corrected for the orbital effects on phase 
residuals, we then showed how this source joins SAX 
J1808.4-3658 (Burderi et al. 2006), IGR J00291+5934 
(Falanga et al. 2005; Burderi et al. 2007) and XTE J1807- 
294 (Riggio et al. 2007) in the class of those AMSPs that 
exhibit a spin up in the 10 12 - 10 13 Hz/s range, with an 
average rate of (3.7 ± 1.0) x 10~ 13 Hz/s. 

We applied different accretion torque models to the ob- 
served phase evolution, but we did not succeed in having a 
significant improvement of its description with respect to a 
constant spin up model either when implying a dependence 
of the spin frequency derivative on the instantaneous accre- 
tion rate, nor when discriminating between different scenar- 
ios of interaction between the accretion disc and the rotat- 
ing magnetosphere. Applying an Alfvenic toque model we 
derived u = (5.6 ± 1.2) x 10~ 13 Hz/s for the spin frequency 
derivative one day after the first observation available. 

The measured value of the spin frequency derivative im- 
plies a peak accretion rate of at least 15% of the Eddington 
limit. Such a high accretion rate indicates that XTE J1751- 
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305, as already noted for IGR J00291+5934 (Burderi et al. 
2007), accretes matter during outbursts at a much higher 
rate than that usually considered typical (a few per cent of 
L Edd ) for AMSPs. 

The equality between the mass accretion rate deduced 
from timing when a simple Alfvenic torque model is consid- 
ered, and the one obtained from a spectral analysis of the 
X-ray emission would place the source slightly farther than 
the Galactic Centre (d ^ 9.1 kpc), although timing-based 
determinations of the distance may be affected by the large 
uncertainties on the assumed torque model. 

We have also discussed the possibility that the mass 
accretion rate is a factor between 1 and 2 higher than the 
one inferred by the X-ray luminosity, if the emission of the 
antipodal polar cap is not visible due to the occultation by 
a thick absorber and re-emission of this energy is out of the 
considered band. An occultation of the secondary cap emis- 
sion is highly improbable to be due to the star itself, while 
it seems more probable to owe to the accretion disc, as GP 
stated that the absence of a significant Compton reflection 
in the emitted spectrum strongly indicates i > 60° . Defining 
the parameter r\ as the fraction of the accretion luminosity 
effectively emitted by the NS in the 0.7 — 200 keV energy 
band, with respect to the one observed (rj — 1 -j- 2), we get 
a distance estimate < 8.5 kpc when r) £ 1.2, while a lower 
limit of ~ 6.7 kpc is obtained in the case rj = 2. 
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